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EQUATIONS OF MOTION OF A ROCKET* 

By F. R. Gantmacher and L. M. Levin 

The equations of motion of a rocket are given in general form; 
a rocket is defined as an apparatus with a liquid or powder rocket 
motor. 


1. THEOREM OF MOMENTUM AND THEOREM OF MOMENTS 

In considering the motion of the rocket, at each instant of time 
only the state of those material particles which at that instant are 
within the control surface passing through the exterior surface of 
the body of the rocket and the exit section of the nozzle shall be 
included. 

In order to obtain the equations of motion of the rocket, the 
following procedure is used. An arbitrary but fixed instant of time 
is considered. A fictitious solid body is denoted by S with 
mass m, which would be obtained if the rocket at the instant t 
solidified and ceased giving off particles . , The solid body S will 
not be homogeneous; in some of its parts, it will have the density 
of a metal and in other parts the density of a gas, and so forth. 

It shall be assumed that the fictitious solid body S is invari- 
ably fixed to the body of the rocket and from the instant t onwards 
(instant of solidification) moves together with the rocket. The 
momentum of the body S shall be denoted by Q. 

The system £ consisting of all the material particles that 
at the instant t entered the composition of the rooket shall also 
be considered. At the instant t the system £ ooincidee with the 
rocket, but at the succeeding instants certain of the particles of 
the system £ will be outside the rocket. The system £ and the 
solid body S have a constant mass equal to the mass of the rocket 
at the instant of time t. The momentum of the system £ relative 


*"0b Uravneniakh Dvizhenia Rakety." Prikladnaya Matematika i 
Mekhanika, Vol. XI, No. 3, 1947, pp. 301 - 312. 
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to an immovable (more accurately Galilean) system of coordinate axes 
will be denoted by K. Let F = dK/dt be the principal vector of 
the external forces acting at the instant t on the rocket (and 
therefore on the system Z ). 

Compare dK/dt and dQ/dt. The motion of each of the particles 
of the system Z shall be considered as compounded. A particle 
moves relative to S (that is, relative to the body of the rocket) 
but the solid body S is transported as a whole. For the absolute 
relative and transport velocities, the notation v a , v r , and v e 
shall be used. Similarly for the accelerations, the notation w a , 
v r , and v e is used. The Coriolis acceleration of the particle 
shall be denoted by J. 

It is noted that v r is equal to zero for the particles of 
the bodyl and for the particles of the powder, whereas w r is equal 
to zero for the particles of the body and those particles of the 
powder, which at the given instant, do not lie on the combustion 
surface. Then^ 

f'XXa =ff - J U-l) 

where J is the principal vector of the Coriolis forces 

and 



The relative velocity of the particle at the instant 
t-j_ = t + dt is denoted by v^. For the elementary change in 


1 In the body are included all fixed particles of the rocket. 

2 In -the following discussion, the subscript i to denote magni- 
tudes relative to the particles will, for simplification, be ommitted; 
for example, v ai and m^ will be written as v a and m, res- 
pectively. 
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velocity Sv r (relative to the body of the rocket)^ , 

5v r = V]_j, - v r = w r dt, hence 

= Z mT lr ' l™ v r (1.2) 

The momentum relative to the body of the rocket at the instant t 
of the particles of the gas at that instant in the rocket is denoted 
by K r and the momentum relative to the body of the rocket of the 
particles of the gas 4 in the rocket at instant t by Kp r . Now 

SK? = K x r - KF 
2nv r = K r 

y, mvpy = Ki r + krdt 

where k r dt is the momentum (in the relative motion) of those par- 
ticles of the gas, which in the time interval dt passed through the 
exit section of the nozzle; k^, is the momentum relative to the body 
of the rocket of the mass of gas passing per second through the 
section of the nozzle or, as used herein, the momentum rate per 
second relative to the body of the rocket; SK?* is the elementary 
change of the relative momentum of the gas occupying a fixed volume 
(within the control surface). 

The vector k^, has the dimensions of a force. 

The force k r is called the equilibrant of the reaction forces 
or simply the reaction force. ® If p = - dm/dt is the mass flow 


> (1.3) 

J 


^Here and in the following discussions, the symbol 5 denotes the 
differential (elementary change) of the vector relative to the body 
of the rocket. The elementary change relative to the initial system 
of axes (fixed) is denoted by d. 

4 Por a rocket with liquid reaction engine, the vector K r is 
the momentum of the particles of gas in the combustion chamber and 
nozzle and the particles of the liquid moving in the tanks and the 
pipes supplying fuel to the combustion chamber. 

^In the equivalent reactive force there are often included 
certain external forces; these will be considered more in detail in 
section 2. 
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per second and u (the average at the exit section), the velocity of 
flow of the gas relative to the nozzle® k r = jiu. From equations 
(1.2) and CL. 3), there is obtained 



&K? 

"ccr 


(1.4) 


Substitution on the right side of equation (l.l) for y.fflWy 
of the expression from equation (1.4) and F for dK/dt yields 

r r 


dQ 

dt 


= F - k + J - 


5K£ 

dr 


(1.5) 


This equation expresses the momentum theorem for the solidified 
rocket, that is, for the solid body S. 

The kinetic moments of the system Z and the body S are now 
considered. 

The notation A, the kinetic moment of the system Z , 
and L of the body S in the absolute motion, that is, in the 
motion relative to the fixed system of the coordinate axes, are 
introduced. The pole relative to which the kinetic moment is taken 
is denoted by a subscript; thus, for example, L c is the kinetic 
moment of the body S relative to its center of inertia C and A q ^ 
is the kinetic moment of the system relative to its center of 
inertia Cq. 

Together with the absolute motion, the motion relative to axes 
passing through the point C and moving forward together with it 
must be considered. 6 7 The magnitudes referring to this motion shall 
be denoted by a prime. 


Similarly in considering the motion relative to axes passing 
through Cq (the center of inertia of the system Z ) and having 
a forward motion the corresponding magnitudes will be denoted by 
a double prime. 


6 

If the rocket has several nozzles, k r = /muq, where |iq 
is the gas flow of the i th nozzle and uq is the mean velocity 
at the exit section of this nozzle. In the following discussion, 
a rocket with a single nozzle is considered; this assumption does 
not affect the generality of the results. 

7 The origin of this system of coordinate axes C is not dis- 
placed relative to the body of the rocket. 
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Applying the theorem of the moment of momentum relative to the 
center of inertia to the system Z , there is obtained for the 
instant of time t 


fv; 

at 


= Gr, 


Cl 


( 1 . 6 ) 


where is the principal moment of all the external forces acting 

on the rochet (and therefore on the system S ) at the instant of 
time t. 

It is notea that A c " =A c f . For, in passing to another system 

of coordinate axes moving translationally relative to the first, there 
is added to the velocities of all points of the system the same 
velocity constant in magnitude and direction. The additional 
momentums will he proportional to the masses and are in the same 
directions. Hence, they reduce to a single equivalent resultant 
vector applied to the center of inertia C^. The moment of this 
additional momentum vector relative to C-j_ . will he equal to zero. 
Further 


A c i = Ac + C-jO + IC* 


whence 


**01* dA c' + dc l c 


cLK 


dt 


dt 


r xr + c i c x ^ 


dt 


(1.7) 


because at the instant t the point coincides with the 

point C and 

dC-jC i » 

— L- = v c - v cl » - — K 
dt c ci m 


By noting that at the Instant t the points and 

coincide and therefore = G- c and recalling that A C ]_" 

there is obtained from equations (1.6) and (1.7) 

dA 1 

-=£- = G 
dt c 

Again the motion of each particle of the system £ shall be 
‘considered as compounded. The motion of the particle relative to the 
axes moving forward together with the point C shall be considered 
as absolute, the motion of the body S relative to these axes as the 
transport motion, and finally the motion of the particle relative to 


= Aci' 
U-8) 
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the body S (that is, relative to the body of the rocket) as relative 
Then 


A c ’ " I r x 

dA. ’ Ti . 
— - 2 rx 


mv„ 


mw 


where r is the radius vector of the particle drawn from the point 
The sum w r + w e + J is substituted instead of the absolute accel- 
eration of the particle and it is noted that 

dl ’ 


mw~ 

e dt 


Z»x 

- Jr X mj = H c 


(1.9) 


where is the principal moment of the Coriolis forces. There is 

obtained 


dA ' 

c 

dt 


dL c ' 

dt 


- E c + 


Jr x 


row v 


( 1 . 10 ) 


Again let be the relative velocity of the particle at the 

instant t]_ = t + dt. Then Vj^ - v r = w r dt and therefore 


Jr x mwydt = ^r X mv^ - ^r x 

It is noted that 

Jr X mv r = A c r 


32V v 


(l.u) 


( 1 . 12 ) 


where A c r is the kinetic moment of the gas within the rocket in the 
relative motion® at the instant t. The -value of this kinetic 
moment of the gas at the instant of time ti is denoted 
by A lc r • Further 


x mv 


lr 


= Six ravjj. - Jjr ! - r) X v^ 


8 See the previous note relative to K r . The magnitudes A c r 
like K r refers not to the fixed mass but to a fixed volume 
occupied by the gas. 
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where r-^ is the radius vector at the instant t^ of the particle, 
which at the instant t had the radius vector r. 

From the sum ^r^ x mv^, those components in which the end 
of the radius vector extends beyond the limits of the rocket are 

taken out. The sum of these components will he equal to the kinetic 
moment per second of the gas through the exit section of the nozzle 
In the relative motion Z rc multiplied by dt. 

The sum of the remaining components gives A, 

xc 

Further the sum ^( r i - r) X mv^, = 0, because 
r l - r = Sr = v r dt and with an accuracy up to infinites- 
imals v^ ~ v r . Thus 

S' X mv^ = A lc r + l rc (1.13) 

X 

From equations (l.ll), (1.12), and (1.13), 

&A r 

X mw = — — + l (1.14) 

r at rc 

r r r 

where 5 A q = A-^ c . - A c is the elementary change relative to the 

body. 

By substituting in equation (1.10) in place of the sum its 
expression from equation (1.14), there is obtained 



aA 


c 


I 


dr 


cLL 1 
c 

dt 


&A c r 

+ l rc - H c + - 


(1.15) 


The kinetic moment of the gas per second in the relative motion 
Z r0 has the dimensions of the 'moment of a force. The moment Z rc 
is called the reactive moment. From equations (1.8) and (1.15), there 
is obtained 


✓ 



■'re 


+ H„ 


5A c r 

dt 


(1.16) 


Equation (1.16) determines the derivative with respect to time 
of the kinetic moment of the solidified rocket S in its motion 


^In the reactive moment, there are often included additional 
moments of certain external forces. See section 2. 
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relative to axes passing through the center of inertia C of the 
forward moving body S. 


2. REACTIVE FORCES 

The area of the exit section of the nozzle is divided into 
elementary areas d(j. The mass per second flowing through the 
area da is denoted by vda and the relative velocity of the 
gas passing through this area by v r . The mass of gas per second 
passing through the exit section of the nozzle will then 
be |i «= Jvdc, where the summation is taken over all the 
elements da of the exit section. 

The vector - tfdcv r , having the dimensions of a force, are 
considered. This force is called the elementary reactive force. 

It originates in the particles of the gas separating from the 
rocket through the area da. The principal vector of the elemen- 
tary reactive forces is equal to - k r , and the principal moment 
relative to the pole C is equal to - Z rc where kj, and l rc 
are the momentums per second and kinetic moment of the gas relative 
to the body of the rocket. 

In the system of reactive forces there are, however, often 
included certain external forces, namely, those arising from the 
atmospheric pressure on the body of the rocket and from the pressure 
of the Issuing parts of the gas on those remaining in the rocket, 
and additional forces due to the unsteady motion of the gas. This 
phenomenon is made clear by the following considerations. 

The combustion in the rocket apparatus is assumed to be at 
standard conditions with the rocket immovably fixed. In this case 

Q = 0 H c = -0 

L c * = 0 F = F* + F* 

J = 0 & 0 - G c * + G c * 

where F # and G #q are the principal vector and the principal 
moment of the forces of the atmospheric pressure and pressure of the 
external part of the gas on the exit section of the nozzle. 
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From equations (1.5) and (1.6), then 10 

SK 17 


- F # - - k r + F* - 


*o 


- i + g * - 

rc c 


at 

sA r 

o 

at 


In measurements on the test stand there are generally deter- 
mined the pressure forces of the rocket on the supports. These 
pressures are characterized hy the principal vector -F and the 
principal moment -G #c , which include in addition to the purely 
reactive forces the additional forces and moments F*, G c *, 
-8K r /dt, -bA r c /dt . 


In the following discussion into one system shall he com- 
bined: (l) the purely reaotive forces, (2) the forces arising 

from the atmospheric pressure 11 and the pressure of the external 
part of the gas (issuing from the rocket jet), and (3) the addi- 
tional forces due to the nonsteadiness of the motion of the gas 
in the rocket. All these forces will be included in the system 
of reactive forces. The principal vector of these forces -F* 
is denoted by T and the principal moment -G* 0 by Mo* Then 


T = - k r + I* - 

sir 

M c - “ *rc + °c* " 


( 2 . 1 ) 


Generally in computing the reactive force T and the reactive 
moment M q the third components in equations (2.1), that is, 

&K r /dt and &A r 0 /dt are neglected. 12 


°It is here assumed that the flow of the gas in the chamber and 
nozzle for the moving and stationary rocket is the same. This 
assumption is equivalent to neglecting the effect of the acceleration 
of the rocket on the relative motion of the gas* 

11 By the forces of the atmospheric pressure is meant the forces 
due to constant atmospheric pressure on the external surface of the 
rocket at standard conditions. 

12 That is, the motion of the gas in the rocket is considered as 
quasi-stationary. 
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3. FINAL FORMULATION OF THE FUITOAMENTAL THEOREMS. PRINCIPLE OF 

SOLIDIFICATION 


From the number of external forces were separated the forces 
arising from the uniform atmospheric pressure and the pressure of 
the external gas, which were included (section 2) in the system of 
reactive forces. These forces have the principal vector F* and 
the principal moment G c *. The principal vector and the principal 
moment of the remaining external forces are now denoted by F 
and G. 


With this notation, the equation of the change of momentum (1.5) 
and the equation of the kinetic moment for the solidified rocket (1.16) 
assume the form 


dQ 

dt 


= F + F* - 


v t 

*r + J - dT 


aV 


dt 


G c + G c* 


ler + H c 



By introducing these equations the reactive force T and the 

moment M according to equations (2.1), there is finally obtained 
c 


dQ 

dt 


<= F + T + J 


dL • 
c 

dt 


= G + M + H„ 
c C G 


(3.1) 


From these equations the following principle of solidification is 
obtained: 


The equations of motion of a rocket at an arbitrary instant of 
time t may be written in the form of the equations of motion of 
a solid body of constant mass if it is assumed that the rocket became 
rigid and solidified at the instant of time t (that is, ceased to 
give out particles) and that to the fictitious solid body thus 
obtained there were applied: (1) external forces acting on the 

rocket, (2) reactive forces, and (3) Coriolis forces. 
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4. SECOND DERIVATION OP -THE EQUATION OF MOMENTUM AND THE EQUATION 

OP MOMENTS 

An Infinitely small interval of time 4t is considered and the 
elementary increments, for this interval, of dK and dQ the momen- 
tum of the system £ , and the solid hody S are compared. 

For this purpose, the volume oocupled by the rocket is 
divided into three parts: (1) the volume occupied by the body, 

(2) the volume occupied at the instant t by the powder (or liquid) 
fuel, and (5) the volume occupied at the instant t by the particles 
of gas in the rocket. The parts of the momentum referring to these 
volumes are denoted by the subscripts 1, 2, and 3, respectively. 

Evidently for the body, 


^ " 4Q 1 (4.1) 

The powder is next considered. At the instant t, K = Q^j 
but at the instant t 1 = t + dt the vector K 2 differs from Q 2 

by the momentum of the elementary mass of powder that burned in the 
time interval dt. This deficiency of momentum is presented in the 
form kpdt where k^ is the momentum expenditure per second of 
the powder. Then^ 

dEg = dQg - k p dt (4-2) 

The volume of the part of the rocket occupied at the instant t 
by the gas is now considered. The momentum of the particles of gas 


^ 3 The momentum of the gas occupying at the instant t^ the volume 

of burned powder after dt seconds is not considered here. The 
ratio of this momentum to the value kpdt is equal to the ratio of 

the density of the gas to the density of the powder and is very small. 
Thus by k_, the expenditure per second of the momentum of the powder 
multiplied P by 1 - e is known. For the rocket with liquid 
fuel dKg = dQg - k^dt + dK^ r where K^ r is the momentum of the 

liquid fuel relative to the motion and Is: the consumption 

per second of the momentum of the fuel irr the transport motion mul- 
tiplied by 1 - € , where € is the ratio of the density of the 
fuel vapors to the density of the fuel itself. 
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in this volume relative to the "body of the rocket is denoted (as in 
section 1) by K r . The corresponding increment is dK r , vhere the 
differential d takes account of the vector relative to the sta- 
tionary system of coordinates. But then 

dK 1 * = SIC 1 * + 4 ) X K 1 dt 

where 6K r is the elementary change of the vector K r relative to 
the body of the rocket and Ci> is the angular velocity of the body. 


The increment of momentum of the gas in the volume considered 
in the transport motion is now considered. This increment agrees 
with dQj in the case where at each point -of the volume considered 
the density of the gas p at the instants t and t^_ is the same. 
In the general case, however, it is necessary to take into account 
also the change in density p. Hence, in the general case the 
required increment in the transport motion is equal to 


4Ci 3 + W VT 


where the integral is extended over that part of the volume V that 
at the instant t is occupied by the gas, v e is the transport 

velocity of the element of volume dT, and Sp/5t is the rate of 
change of the density p in the given volume element dT. Hence 


= dQ 3 + 



St 


dtv e dT + 5K r + (i) x if dt 


(4.3) 


The increments in the momentum of the particles within the 
volume of the rocket have been considered. In computing dEC, 
however, it is necessary to take into account also the moment at the 
time t-^ of those particles, which in the interval from t 

•to t-j_ = t + dt passed through the exit section of the nozzle. 

This momentum is equal to (k e + kp)dt where k e is the rate of 
expenditure per second of the momentum of the gas through the exit 
section of the nozzle in the transport motion and kp in the 
relative motion. Thus 


+ dKg + aitj + (k e + kp)dt 
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Substituting dKL^, dEC,, and cLK^ from equations (4.1), (4.2) 
and (4.3), respectively, there is obtained 


dK 


= dQ - (kp - k e - kj.)dt + 61^ ttoxf dt + |£ dt v dT (4.4) 


As before, the principal vector of the forces of the atmospheric 
pressure and the pressure of the external part of the gas on the 
rocket is denoted by F* and the principal vector of the remaining 
external forces by F. According to the momentum theorem for the 
system I! , cLK/dt = F + F* and therefore from equation (4.4) 
there is obtained 


dQ 

dt 


= F 




k 

P 


k 

e 



v dT - <0 x K* 
3t e 


(4.5) 


where the last term - x K r represents half the Coriolis force J. 
Hence, by taking equations (2,1) into account, equation (4,5) can he 
reduced to the form 

= F + T + k -k - \\\ dT + - J (4.6) 
dt P e \\\ St 6 2 


Similar considerations permit obtaining the equation of moments 


dL’ 

ST " G o +M ° + W xA ° r 


r' x v e ' 



(4.7) 


where Z c p is the kinetic moment of the powder per second and 2 ce 
is the kinetic moment of the gas per second through the exit section 
of the nozzle. In both cases, there is considered the kinetic moment 
in the transport motion of the body relative to a system of axes 
moving forward together with the center of inertia C of the body S. 


•*- 4 This equation holds also for the liquid fuel rocket. In this 
case, however, K r is the momentum of the gas and the fuel in their 
relative motion. 
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5. EQUATION FOR THE MOMENTUM, CORIOLIS FORCE, AND CORIOLIS MOMENT 

By comparing equations (4.6) and (4.7) with equations (3.1), 
the equations for the Coriolis force J and the Coriolis moment H c 
are obtained: 



These equations are evidently simplified if the density of the gas 
does not vary with time, that is, dp /at = 0. (This equation holds for 
the quasi- stationary motion of the gas in the rocket.) 


The expressions (5.1) for J are first considered. The mass of 
powder expended per second (or what amounts to the same thing, the 
addition of the mass of gas) is denoted by p^ and the mass of gas 
per second through* the exit section of the nozzle by p 2 » Then 
evidently 


“l " 


V"? ~ ^ 4 - 



p dT = 



dT 

St 


(5.3) 


Further, k^dt represents the momentum of an elementary mass 
of powder p-^dt burning in the interval of time dt. The center of 


inertia of this elementary mass of powder is denoted by C^ and 
the velocity of the point Cp in the transport motion of the body 
by v cp . Then kp = PjV^. 


If it is assumed that the elementary mass of powder p^dt is 

symmetrically placed relative to the axis of the rocket, the point 

C will lie on this axis. If the mass p,dt is symmetrical with 
P 1 

respect to the mean section of the powder containers, C will be in 
the plane of this section. ^ 

In the same manner, the elementary mass of gas p^dt passing 
through the exit section of the nozzle in time dt is considered. 

By C q the center of inertia of this elementary mass is denoted and 
by v ce the transport velocity of this center. Then k g = ^ 2 v ce‘ 
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If the gaa passing through the exit section Is symmetrical, the 
point C e lies on the center of the exit section of the nozzle. 

By Tp and. r Q the radii vectors of the points Cp and C Q , 
respectively, drawn from any pole 0 fixed to the body are denoted. 


Let r denote the radius vector of an arbitrary point within 
the rocket. By using the expression for the velocities of the points 
of a solid body. 


▼op - T o + W x r 
T c0 = v 0 +0) x r 
v e = v Q + <0 x r ^ 



(5.4) 


On the "basis of these equations and the equations k^ = ^q v cp 

and k e = \i^ v ce> equation (5,1), taking account of equation (5.3), 
assumes the form 


J 





(5.5) 


By comparing this equation with the usual expression for the 
Coriolis force J = - 2 <*) X K r and taking into account the arbi- 
trariness of the vector 0), there results 


K r = U 2 r e - Ui r p + 


3p 

Bt 


tt- r dt 


(5.6) 


where K r is the momentum of the gas in the rocket in its motion 
relative to the rocket and the integral on the right side is 
extended over the entire volume W occupied by the gas. 


Equation (5.6) is of a general character.-^ It determines the 
momentum of the gas (or liquid) enclosed in the given volume in the 


15 Equations (5.6) and (5.2) for X r and H c can be obtained 
directly if the derivatives are computed with respect to. time of the 
integrals 


t. 



rdm r x (a) X r) dm 

taken over the mass occupying the given volume at the time 
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case where the incoming and outgoing masses p-, and pg of the gas 
through the surface hounding this volume and the per second change in 
density at each point of the volume are given. In this equation, rp 
is the radius vector of the center of inertia of the incoming mass per 
second pp and r e the radius vector of the outgoing mass pg. 

The particular case where dp/dt 2 0 (for example the quasi- 
stationary motion of a gas, the arbitrary motion of an incompressible 
liquid-*-®, and so forth) is considered. In this case according to equa- 
tion (5.3), there is pp = pg = p where for the case of the rocket 
p = dm/dt. If the vector b = C e Cp is introduced, equation (5.6) 
and the equation for J become 

K r =pb J = -2(0xK r = -2pfo>Xb) (5.7) 

Equation (5.2) for the Coriolis moment H c is now considered. 

If it is assumed that 5p/St = 0 and that the equivalent vector of 
the momentum of the particles of gas in their relative motion passes 
through the center of inertia C, that is, A Q r = 0, equation (5.2) 
"becomes 


H = 7 1 - 7 » 

^c 1 pc L ec 

Eor example, the case of plane-parallel motion of the rocket is 
considered. The moment of inertia of the body relative to the equa- 
torial axis passing through the point C and perpendicular to the 
plane of the motion is denoted by I. The moment of inertia of the 
mass of powder burned in dt seconds will then be equal to -dl. 
Hence the kinetic moment of this mass will be 

! 'po dt - - ““ 

On the other hand, the kinetic moment (in the translational 
motion) of an elementary mass of the gas passing through the exit 
section in time dt will be -dm r 2 (i) where r 0 = UU 0 . Thus 
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6. EQUATION OF MOTION OF THE CENTER OF INERTIA OF THE ROCKET 

From the principle of solidification for the instant of time t, 
there results 


xnw c = F + T + J (6.1) 

where m is the mass of the rocket at instant t and v c is the 
acceleration of the center of inertia C of the solidified rocket 
at time t, that is, of the "body S. 

During the combustion of the fuel, however, the center of 
inertia of the rocket is displaced relative to the body. The motion 
of the center of inertia of the rocket relative to the initial 
(fixed) system of coordinates can be represented in the form of a 
compound motion in which the center of inertia moves relative to the 
body (relative motion) and the body of the rocket moves relative to 
the fixed system of axes (transport motion). Then v c and w c 
are the transport velocity and acceleration and v or and w cr the 
relative values. The absolute velocity and the acceleration of the 
center of inertia are determined from the equations 

v = v 0 + v cr 


w = w c + w or + 2(0 X v cr 


Determining from the second equation w 0 and substituting in 
equation (6.1), there is obtained the equation of motion of the center 
of inertia of the rocket. 

mw=F+T + J + mw cr + 2m d) X v cr (6.2) 


The expressions for the magnitudes v cr and w cr shall be 
found. The body of the rocket may here be considered as fixed. Let 
the point 0 of the body of the rocket be taken as the initial 
and C and 0* the position of the center of inertia of the rocket 
at the instant t and t * = t + dt . Then for the moments of the 
times t and t’, there results 


mr 


c - 


(m + dm) r c ' 



p (t + dt) r dT 
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where r c = W, r c i = TXJ* = r c + 5r c , and p is the density of the 
elementary volume of the rocket. The integration is extended over 
the entire volume of the rocket (including the hody). 


By subtracting the first relation from the second and neglecting 
the terms of the second-order smallness, 


dm r c + m 5r c = - p-]_ r p + 




dt r dT 

dt 


’whence it follows that 


5r 


IQY, 


cr 


= m 1? ° ^2 r c “ 1*1 r p 


r- + 



dp 

St 


r dT 


or on the basis of equation (5.6) 

m ^cr = ^2 ( r c “ r e ) + ^ 
By differentiating the equation (6.3), 

dp. 


SK 


or 


- ^2 v cr + ^cr = ^2 T cr + -§£ (r c “ r e> + Kt 
mw cr = 2p 2 v or + — (r c - r 0 ) + 


dt 


(6.3) 


(6.4) 


In the quasi-stationary case dp /St CO, Pg = p-^ =• p, 

5K r /dt » 0, dpg/dt *' 0, and equations (6.3) and (6.4) assume the 
form 


v H. (r - r ) = ii c G 

cr m ' o P' m P 


(6.5) 


w cr = 2 ^ ( r c - r p) = 2 S ^ 

m id c * 


The computations conducted by these equations show that the mag- 
nitudes v, and W-™ are negligibly small in comparison with the 

CjT vr 

mean velocities and the accelerations on the active part of the 
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trajectory. Hence in equation (6,2), the term nrw cr may "be neglected 
in comparison with the reactive thrust T. The term 2 m a) X v or is 
generally of the same order as J because | r c - r ? | is of the same 
order as |r e - rp |. 

7 . EQUATIONS OE THE ROTATIONAL MOTION OF THE ROCKET 

By £, r) , and £ are denoted principal central axes of 
inertia of the rocket, by Ip, Ig, and I 3 the' moments of 
inertia of the rocket relative to these axes, by p, q, and r 
the projections of the angular velocity of the body ( 1 ) on these 
axes. 


It is first assumed that the directions of the axes of iner- 
tia | , 7 ], and £ are fixed relative to the body during the 

entire time of combustion. This natural assumption is always made 
for rockets because it is assumed that at all times the axis of the 
rocket is one of the principal axes of inertia and the two other 
principal axes of inertia may be arbitrarily chosen in the plane 
perpendicular to the axi3 of the rocket. 

Under the assumption made, the principal axes of inertia of the 
solidified rocket, that is, of the body S, will at all times be 
parallel to the axes | , q, and £. The corresponding moments of 
inertia for the body S will have constant values equal to the 
values of the moments of inertia Ip, Ig, and I 3 of the rocket 
at the instant of time t. It follows that p, q, and r are the 
projections of the angular velocities of the solid body S on the 
principal central axes of inertia of this body. Hence, making use 
of the principle of solidification, the three equations of the rota- 
tional motion can be written in the form of Euler 17 

II —E + (I 3 - ± 2 ) qx* = N| (| , H, £ > 1* 2, 3; p, q, r) 

(7.1) 


17 Here and in the following discussion, the parentheses after 
each equation indicate that two other equations are obtained by 
cyclical interchange of the letters and subscripts in the paren- 
theses. 
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where Ng is the sum of the moments relative to the axis £ of all 

the external reactive and Coriolis forces and similarly for 
and Ng . 

The three equations (7.1) together with the three scalar equa- 
tions of motion of the center of inertia constitute a system of six 
equations, which determines the motion of the rocket. 


For completeness, the general case when the principal central 
axes of inertia change their directions relative to the body of the 
rocket shall he considered. 

By £ 1 , t) ' , and £ ' the directions of the principal axes 
of inertia of the body S is denoted, that is, the directions fixed 
relative to the body and coinciding with the directions £ , q, 
and £ at the time t. The projections of the angular velocity to 
on the axes £ ', q‘, and £* are denoted by p 1 , q', and r* 

respectively. Applying the principle of solidification, for the 
moment of time t 

II Jr + (I 3 - *2> q' r* = ( £, T|, £ j 1, 2, 3; p, q, r) 

(7.2) 

where Ng' is the sum of the projection of the external reactive 

and Coriolis forces on the axis £', and so forth. By £2 the 
angular velocity of the trihedron £q£ relative to the body of 
the rocket is denoted. Let dto/dt be the derivative of the vector 
relative to the body of the rocket and Sto/dt the derivative 
relative to the axes £q£ . Then 

~ = !?+ ft X to (7.3) 

dt dt 

The projections of the vector dto/dt on the axes £’, q’, 
and £ 1 are equal to dp’/dt, dq'/dt, and dr'/dt. The pro- 
jections of the vector 5to/dt on the axes £ , q, and £ are 
equal to dp/dt, dq/dt, and dr/dt. Because at the instant t 
the directions £ 1 , q ' , and £ 1 coincide with the directions 
q, and £ for this instant of time, from equation (7.3) 

tt“ = It + " ^£ g ( t > n> Py q> r ) 


(7.4) 
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It Is noted that p’ = p, q' = q, and r ' = r at the 
instant t. Substituting in equation (7.2) in place of dp ' /dt , 
dq ' /dt , and dr’/dt their expressions from equation (7.4) 
there is finally obtained 


5 


x i If + (I 3 - J 2 ) <i r + J i <Pn r " = w e 


(| y "H* t *> Pj 


(7.5) 


where fig, and Q* must "be considered as known functions 

of the time. * 
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